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Abstract 

It is shown that the exponential moments of the canonical phase can be directly 
sampled from the data recorded in balanced homodyne detection. Analytical 
expressions for the sampling functions are derived, which are valid for arbitrary 
states and bridge the gap between quantum and classical phase. The recon- 
struction of the canonical phase distribution from the experimentally determined 
exponential moments is discussed. 



1 Introduction 



Since Dirac's attempt in 1927 to introduce amplitude and phase operators in 
quantum mechanics a number of concepts have been developed with the aim 
to overcome the problems resulting from the non-existence of a Hermitian phase 
operator (for a review, see Recently an attempt has been made to bridge the 
gap between two concepts which are based on essentially different approaches to 
the phase problem and widely used in quantum optics P]. In the first, the phase 
of a radiation-field mode is defined from the requirement that phase and photon 
number should be complementary quantities. This first-principle definition leads 
to the canonical phase (also called London phase), the associated phase states 
being the right-hand eigenstates of a one-sided unitary exponential phase oper- 
ator 1^. In the second, phase quantities are defined from the output observed 
in phase-sensitive measurements, such as eight-port homodyne detection. It is 
well known that in such a scheme the Q function [or, in the case of non-perfect 
detection, a smoothed Q function, i.e., an s-parametrized phase-space function 
with s < —1] is measured P, P). The measured phase distribution can then be 
obtained from radially integrating the (smoothed) Q function. Whereas in the 
classical limit the measured phase coincides with the canonical phase, in the 
quantum regime the two phases significantly differ from each other in general, 
because of the additional noise unavoidably connected with the Q function. So, 
from a study of the asymptotic behaviour of the measured and canonical phase 
distributions in the semiclassical domain it can be anticipated that the measured 
distribution is at least broader than the canonical one . In the quantum regime 
it is principally not possible to obtain the canonical phase distribution from the 
radially integrated Q function, but it must be related to the complete quantum 
state, i.e., the complete Q function as a representation of the state in the phase 
space. 

The best and perhaps ultimate method for measuring the quantum state of a 
traveling optical field has been four-port homodyne detection in which the quan- 
tum state is measured in terms of the quadrature-component distribution ||^. 
Since the quadrature-component distribution contains all knowable information 
on the quantum state, the various quantum-statistical properties of the system 
can be obtained from it. Moreover, the quadrature-component distribution is 
less noisy than the Q function and therefore it is more suitable for determin- 
ing the quantum statistics than the Q function. The method also called optical 
homodyne tomography (OHT) was first used for reconstructing the Wigner func- 
tion of a single-mode optical field applying inverse Radon transform p[, which 
requires a three-fold integration of the measured data. In the numerical calcula- 
tion the standard filtered back projection algorithm is usually used, so that the 
reconstruction of the Wigner function is biased by data filtering. 

This problem does not appear and the effort is drastically reduced if the 
quantities that one is interested in can be directly sampled from the homodyne 
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data. In particular, the determination of the quantities and the error estimation 
are very fast and can be performed in real time. Systematic errors can easily be 
reduced to any desired degree of accuracy and the remaining statistical errors only 
reflect the finite number of measurement events. It has been shown that both the 
density matrix in the photon-number basis and the moments and correlations 
of the photon creation and destruction operators |]10| can be obtained in this 



way, which has offered novel possibilities of the experimental determination of 
the photon-number statistics of light. 

In contrast to the photon-number statistics, the determination of the canoni- 
cal phase statistics has been an open problem. The phase statistics can of course 
be tried to be determined indirectly by calculating it approximately from the 
Wigner function, using in the reconstruction of the Wigner function the standard 



filtered back projection algorithm of inverse Radon transformation [11]. Another 



indirect method, which avoids the rather lengthy detour via the Wigner func- 
tion, is the calculation of the phase statistics from the sampled density matrix 
in the photon- number basis |T^. However, since canonical phase and photon 
number are complementary quantities, there is no a priori upper bound for the 
density-matrix elements that can contribute to the phase statistics. Hence, large 
numbers of density-matrix elements must be sampled, the statistical errors of 
which then give rise to an error accumulation in the phase statistics such that 
the inaccuracies eventually dominate the result (note that the statistical error 
of the off-diagonal density-matrix elements increases with the distance from the 
diagonal). To limit the effect of inaccuracies, one must necessarily restrict the 
method to states of low photon numbers and appropriately truncate them. A way 
that remains to overcome the problem is to directly sample the phase statistics 
from the homodyne data. Unfortunately the canonical phase distribution cannot 
be related to the quadrature-component distribution in the sense of a sampling 
formula because the corresponding integral kernel does not exist. It has been 
therefore suggested to introduce the exact phase distribution as the limit of a 
convergent sequence of appropriately parametrized (smeared) distributions each 
of which can directly be sampled from the homodyne data [1^ . The exact phase 
distribution can then be obtained asymptotically to any degree of accuracy, if the 
sequence parameter is chosen such that smearing is suitably weak. In practice 
it is therefore required that the sampling procedure is performed simultaneously 
for various values of the sequence parameter, each value giving rise to its own 
sampling function. The disadvantage of the method is rather technical, since the 
numerical effort drastically increases with the number of photons contained in 
a state. This fact makes the method effectively applies only to states with low 
photon numbers. Finally, it has been suggested to measure the canonical phase 
distribution asymptotically by replacing the local-oscillator in the homodyne de- 
tection scheme with a reference mode prepared in so-called reciprocal binomial 
states - a method that is also state dependent and hardly realizable at present 

m. 
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In this paper we show that the problem of direct determination of the canon- 
ical phase statistics from the homodyne data can be solved when it is based on 
the exponential phase moments (i.e., the Fourier components of the phase distri- 
bution) and not on the phase distribution itself. We show that the exponential 
phase moments can be directly sampled from the quadrature-component distribu- 
tion, without making a detour via other quantities and without any assumptions 
and approximations with regard to the state. In particular, we derive analytical 
expression for the sampling functions and give a very simple procedure for the 
numerical calculation. Since the method is independent of the state, it applies to 
both quantum and classical fields and all fields in between in a unified way and 
bridges, through the universally valid sampling functions, the gap between quan- 
tum and classical phase. Needless to say that for obtaining the full information 
on the phase statistics, all (non-vanishing) exponential phase moments must be 
determined. It is worth noting that already sampling of a few low-order moments 
provide us with interesting information [jl5| . 



The paper is organized as follows. In Sec. ^ the problem of direct sampling 
of the exponential phase moments of a classical oscillator from the quadrature- 
component distribution is studied. In Sec. § the theory is extended to the canon- 
ical phase of a quantum oscillator. Measurement errors are studied in Sec. ||, 
and in Sec. ^ numerical results of computer simulations of measurements for de- 
termining the canonical phase statistics are presented. Lengthy mathematical 
derivations are given in appendices. 



2 Sampling of exponential phase moments — clas- 
sical case 

In order to gain insight into the problem of phase measurement by means of 
balanced homodyne detection, let us first consider the situation in classical optics. 
Here we can assume a proper phase-space probability W{q,p) dqdp, which can be 
rewritten as, on introducing polar coordinates q = r cos if and p = r sin ip, 

W{q, p) dqdp = P{r, ip) drdip, (1) 

where 

P{r,(f) = rW{r cos (p,r simp). (2) 
The phase probability distribution P{ip) is then defined by 

P{ip)= / drP{r,ip), (3) 
Jo 

and the exponential phase moments \I'fc, which are given by the Fourier compo- 
nents of the phase probability distribution, 

■^k = I dp>e'^^P{p>), (4) 

J2tt 
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can be written as 



= dip dre'^'PP^r.ip). (5) 



27r 



In balanced homodyne detection the probability distributions d) for the field 
quadratures x{d) =qcosd + psmd are measured. [Note that when the harmonic 
oscillator represents a moving particle in a harmonic potential, then d = ujt is 
valid and x{ujt) is the time-dependent position of the particle.] The quadrature- 
component probability distribution d) can be obtained from the phase-space 
probability distribution P{r,ip) as, on recalling that =r cos{ip—{}), 

p{x,d) = dip dr P{r, cp) 5[x — r cos{(p—'&)], (6) 

J27r Jo 

which is nothing but the well-known Radon transform, whose inversion yields 
the phase-space probability distribution in terms of the quadrature-component 
distribution. 

Let us now turn to the problem of direct sampling of the exponential phase 
moments. A quantity A can be determined from the homodyne data by means 
of the sampling method, if it can be related to p{x, -d) as 

A= d§ dxKA{x,^)p{x,^) (7) 

J2tt J —oo 

with a well-behaved integral kernel K_4(x, -(9) as sampling function. Note that 
p{x,{} + vr) = p{—x,{}), so that the {} integration in Eq. (|) can be restricted to 
a TT interval. For the sake of convenience, here and in the following we prefer a 
271 interval. In contrast to the full phase-space probability distribution P{r,ip), 
which cannot be obtained from p{x, i}) by a simple inversion of Eq. in the form 
of Eq. d^), the sampling method applies to the Fourier components of the radi- 
ally integrated phase-space probability distribution, i.e., the exponential phase 
moments can be given by 

f roc 

^k= d§ dxKk{x,'d)p{x,§). (8) 

J2-K J —oo 

Substituting in Eq. (H) for p{x,'d) the expression (H), we arrive at 

/• /■ /"OO /"OO 

^k = j d^ j dip I dx dr Kk{x,^)P{r,ip)5[x-rcos{ip-^)\. (9) 



In order to determine the integral kernel Kkix^d), it is convenient to introduce 
the Fourier decomposition 

oo 

Kk{x,^)= E e^'^^KkA^). (10) 

l=—oo 
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When the phase argument in P{r,ip) is shifted towards ip + ipo, i.e., P{r,ip) 
^P{r,(f + (fo), then it follows from Eq. that changes as \E'fc e~*^'^"\l/fc. 
Comparing this requirement with Eqs. (|) and (0), we find that Kk^i{x) must 
be of the form Kk,i{x) = Sk,iKk{x), and hence 

Kk{x,'d) = e''^Kkix). (11) 

We now insert this expression into Eq. @, compare the result with Eq. (^) and 
find that Kk{x) must satisfy the integral equation 

/ d^e''"^Kk{rcos^) = 1 (12) 

for all r > 0. 

From Eq. (|l^) we can see that Kk{x) is not uniquely defined. First, any 
function of parity (—1)^+^ can be added to K^^x) without changing the integral. 
Second, any polynomial of a degree less than k can also be added to Kk{x). As 
can be verified by direct substitution, a solution of Eq. (|T2|) for odd and even fc, 
respectively, is given by 

K2m+i{x) = |(-l)"^(2m+ l)sign(a;) (13) 

and 

i^2m(a:) = 7r-^(-l)'™+^mln|x| + C (14) 

[m = 0, 1, 2, . . .,where C is an (irrelevant) constant.] Note that this solution 
ensures that the integral (^ exists for any quadrature-component distribution 
p{x,'&) which with increasing |x| decreases at least as |a;|~^"^+'^'', e being a (arbi- 
trarily small) positive constant, i.e., for any physical state. Clearly, this would 
be not the case if, within the ambiguity mentioned, polynomials were added to 
the functions (|13]) and (0). Another reason for choosing the functions (|l^) and 



( p^ is the reduction of the statistical error in a real experiment. Since this error 
is related to the variance of the kernel (Sec. ^I3| ), it is advantageous to choose 
kernels which are varying as slowly as possible. 

3 Sampling of exponential phase moments — quan- 
tum case 

It is worth noting that the results derived in Sec. also remain valid for a quan- 
tized radiation-field mode, provided that W{q,p) [or in polar coordinates, P(r, ip)] 
is identified with the quantum-mechanical Wigner function. Hence, using in 
Eqs. (H) and (|lT]) the functions (|13D and ( p!4D enables one to determine exponen- 
tial phase moments defined by the Fourier components of the radially integrated 
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Wigner function from the homodyne data by means of the samphng method. 
Since the Wigner function of a quantum oscillator can attain negative values, it 
cannot be regarded, in general, as a proper phase-space probability distribution, 
and hence the radially integrated Wigner function does not represent, in general, 
a proper phase distribution function. 

As already mentioned, for a quantized radiation-field mode photon number 
and canonical phase are complementary variables, and in place of Eq. (|]) we have 

P(^) = {2ti)-\^\q\v), (15) 

where g and \ip), respectively, are the density operator of the state and the (un- 
normalizable) phase states |]T6[ 

oo 

|v,) = ^e-1n), (16) 

n=0 

which are right-hand eigenstates of the only one-sided unitary operator 

E={n + lfl^a, (17) 



E\v)=e'^\^). (18) 

In Eq. (p!?!), n = d)a is the photon-number operator, a) and a being the photon 
creation and annihilation operators, respectively. From Eqs. (|15]) and ([T8|) to- 
gether with the fact that the phase states resolve the unity it is easily seen that 
the exponential phase moments defined in Eq. (^ can be written as 

^fc = (^^) (19) 



for /c= 1, 2, . . ., and = ^Ifc for k = —1, —2, . . .. We now combine Eqs. (|T^ and 
( [IQD and obtain [in place of Eq. 



oo 



^k = J2 Qr^.n (20) 
n=0 

(fc = 1, 2, . . .). Next, we express the quadrature-component distribution 

p{x,d) = {x,d\Q\x,d) (21) 

in terms of the density-matrix elements in the photon-number basis. For this 
purpose we expand in Eq. (^Tj) the eigenstates |x, d) of the quadrature-component 
operator x{'d) = 2^^/^ (e~*''a + e*''a^) in the photon-number basis [0, 



oo 

Awd 



a;,^?) = ^e^"V„(x)|n), (22) 



n=0 
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where the functions ipn{x) are the eigenfunctions of the harmonic-oscillator Hamil- 
tonian, ipnix) = {2'^n\y/TT)^^^'^ exp(— x^/2)H„(x), H„(a;) being the Hermite poly- 
nomial. From Eqs. (^) and ( p2D we then obtain [in place of Eq. 



p{x, ^9) = E E ^n(x)V^„.(x)f?^,„e*("-'")'^. (23) 

Let us again assume that can be obtained from p{x, {}) according to Eq. (§). 
Substituting in Eq. (|^) for p{x, {}) the quantum-mechanical expression (p3|) and 
comparing the result with Eq. (|20|) , we see that the kernel Kk{x^d) must be of 
the form ([l^), but now Kk{x) must satisfy the integral equation 

/•oo 

271 dx Kk{x)iJni^,{x)iJn{x) = 1 (24) 



(n = 0, 1, 2, . . .). Equation (p^) plays the same role for a quantum oscillator as 
Eq. (|l^) for a classical oscillator does. From Eq. (p^ and the properties of the 
Hermite polynomials the same ambiguity in the determination of Kk{x) as in the 
classical case (Sec. H) is found. Provided that a Kk{x) exists, any function of 
parity (— 1)'^"^^ and/or any polynomial of a degree less than k can be added to 
Kk{x) in order to again obtain a solution of Eq. (|^). 

We now turn to the problem of construction of an integral kernel that satisfies 
Eq. (p^. For this purpose we return to Eq. (|19D and bring the operator E'' into 
the normally ordered form, 

^ at"exp(-ata)a"+'= 
n=o Jn\{n + k)\ 

oo CO 1 / 1 \m 

= y y ^^^^n+m^n+m+k ^25) 

^0^0 ^n\{n + k)\ 

(the notation : : is used to indicate normal ordering). From the expansion (|25|) to- 
gether with the sampling formula for normally ordered moments and correlations 



of the photon creation and annihilation operators |T0 



27rV2^ 



-1 

^oo 



2-7T 



d'd rfxe^("-'")''H„+„(x)p(x,^?), (26) 



we find after some calculation (see Appendix A| ) that "^k, Eq- QT^), can be written 
in the form of Eq. (H) , 

/" /"OO _ 

^fc=/ d^ dxKk{x)e'^^p{x,d). (27) 

In Eq. (^TD, the integral kernel Kk{x) can be decomposed into two parts, 

Kkix)=Kkix)-Fkix), (28) 
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where for odd and even k, respectively, Kk{x) reads as 

2x(m+l)! 



K2^+i{x) = (-1)' 



and 



K2m{x) = {-ly 



ml 



(2vr)' 



(27r)-+3/2 7o 1 

r2"^$[m + 2,3/2,-x2tanh(rV2)] 



r 



(29) 



m+1 



„2m- 



igmrV2 r$[^+i^ 1/2, -x^ tanh(rV2)] 



sinh'"(rV2) 
and Ffc(x) is the polynomial 



(30) 



Fk{x) 




Hfc_2n(a;) 



i=o V ^ / ^(/ + !)...(/ + A;). 



(31) 



In Eqs. (p9|) and (0), ^{a,b,y) is the confluent hypergeometric function and 
fi'^'^)(z) defined in Eq. (|A 18|) [together with Eq. (|A 17| )j in lAppendi:!^^ can be 
given by power-series expansion ( [Appendix B| ), 



(32) 



m=0 



where 



27r*^/2 



ml 



T{k/2 + m) dx' 



na 



JX) 



(33) 



x=0 



From Eqs. (H) - (|3T|) [together with Eq. (|Fl8|) ] it is seen that Kkix) exists, 
and it can be proved by direct substitution that Kk{x) satisfies Eq. (^). Hence we 
have found a solution of Eq. (^4]) even when the assumption made for constructing 
it fails (i.e., when the moments and correlations (^) do not exist for all m and 
n, see the derivation in [Appendix A| ). It is worth noting that both Kk{x) and 
Kk(x) are solutions of Eq. (plf), because the polynomial Fk{x) in Eq. ( |28| ) reflects 
the above mentioned ambiguity in the solution of Eq. (|^) and can therefore be 
omitted. Further, it can be shown ( [Appendix (J| ) that with increasing the 
solution Kk{x) approaches the classical one, i.e., the asymptotic behavior for 
large |a;| of Eqs. (B^) and (pOP, respectively, is exactly given by Eq. ([T3) and 
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Eq. (|T^. We see that Kk{x) can be used for determining exponential phase 
moments from the homodyne data for all states whose quadrature-component 
distributions p{x,d) asymptotically decrease at least as e > 0, i.e., for 

any physical state. Since Kk{x) applies to both quantum and classical systems 
in a unified way, the sampling method bridges the gap between quantum and 
classical phase. Examples of Kk{x) for various k are shown in Fig. |I|. It is seen 
that they are well-behaved functions, which rapidly approach the classical limit 
and differ from it only in the small region of vacuum fluctuations. 

Let us comment on the numerical calculation of Kk{x) which can be performed 
in a straightforward manner. In particular, from Eqs. (pQ]) and (^), respectively, 
it can be easily proved that ||T5|, [TT 



KAx) 



dt 



<l>(2,|,-x2tanht) 



/o \/i cosh^t 

and [after calculating f2(^^(r^) according to Eq. ( [A 18| )] [ p^ 



(34) 



K2(x) 



27r Jo 



dtlo{t) 



'2t 



sinh t 



1 , 

2 "^(2, 1, -x^tanht 

cosh tsinht ^ 



(35) 



[lo(t), modified Bessel function]. The one-dimensional integrals can then be cal- 
culated numerically using standard methods. In order to calculate Kk{x) for 
arbitrary k it may be convenient to start from Kk{x) as given in Eq. (|A 2|) in 
[Appendix A| and approximate it in a small (nonclassical) interval by a finite sum 
such that 



K,(x) 



(2 



(x) + Fk(x) 



I X I <C Xq^ 



(36) 



and use the classical limit, Eqs. ( pTSD and ([T^), elsewhere. For example, for the 
calculation of the kernels plotted in Fig. |l| it is sufficient to chose the parameters 
/o = 40 and Xq = 4, and to truncate the infinite sum in Eq. ( pT] ) for Ffc(a;) at / = 
103. 



4 Nonperfect detection and measurement errors 

In practice there is always a number of experimental inaccuracies that limit the 
precision with which the exponential phase moments can be determined. In 
this section we restrict attention to three kinds of inaccuracies: data smearing, 
discretization of the phase parameter i), and finite number of measurement events 
(i.e., discretization of a;). 
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4.1 Data smearing 

Since in a realistic experiment the quadrature components cannot be measured 
with infinite precision, we may assume that instead of p{x, a smeared distri- 
bution 

/oo 
dy f{x -y;T])p{y,i9) (37) 
-oo 

is obtained. In Eq. (pT)), f{x;ri) is some positive single-peaked function and rj is 
a parameter quantifying the smearing effect. A typical example of /(x; rj) is a 
Gaussian, such as 

-7ip-p(-|l). .= ^1^. (38) 

which corresponds to the use of nonperfect photodetectors whose efficiency t] is 
less than unity. 

Substituting in Eq. (||) for the exact distribution p{x, -d) the smeared distri- 
bution p[x,i)]rj) yields exponential phase moments "^kiv) that differ from in 
a systematic error (bias) A^^^''^^ as follows: 

^k{Tl)=l d§ r dxKk{x,^)p{x,§-Ti) = mk + A^^^k (39) 

J2'iT J —oo 



with 



where 



A(^)*fc= / dde'^^ dxgk{x;r])p{x,^), (40) 

J2n J — oo 



9k{x;r])= dyKk{y)[f{y-x;r]) ~6{y-x)]. (41) 

J — oo 

Examples of the kernel gk{x;ri) for the determination of the systematic error 
A^^^^fc are plotted in Fig. |^ [with f(x;ri) according to Eq. (|38D]. From a com- 
parison of gk{x; rf) with Kk{x) (see Figs. |I] and ^ it is expected that the absolute 
values of '^ki'n) smaller than those of ^k in general. The systematic error 
is state-dependent as it can be seen from Eq. (^0|). To give an impression of 
its magnitude, let us restrict attention to the classical limit and consider a state 
whose phase-space probability distribution is radially sharply localized at r = 
ro such that ro ^ a. In this case it can be shown that, on assuming Gaussian 
smearing and using the results in Sec. ||, 

^'fc(//) «iexp( -^j ^fc, (42) 

which reveals that the exponential phase moments can be determined from the 
smeared data quite rehably as long as k<^ro/a. 
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It is worth noting that under certain circumstances it is possible to com- 
pensate for the systematic error during the samphng process, introducing an 
appropriately modified kernel Kk{x\rj). Let us again assume Gaussian smearing, 
which is typically observed in nonperfect detection, and apply Eq. (|37|) together 
with Eq. (pHI). In this case we may replace Eq. ( p6D with ||18[ 

|- / I -| ~i 



\ra+r7i 



n+nv 



X dxe'^"^^}in+m{x)pix,'d;ri), 

J2tt J—oo 



(43) 



and follow the lines given in Sec. |^ and [Appendix A] . It is easily seen that in 
Eq. (p7|) p{x,{}) and Kk{x), respectively, must be replaced with p{x,{};ri) and 
Kk{x;r]), provided that Kk{x;r]) exists. The kernel Kk{x;ri) obviously compen- 
sates for the losses associated with nonperfect detection and can be obtained 
from Eq. (1^1 ) in [Appendix A| , if c/^^ is replaced with ^(r^) = r7-(^+'=/2)cp). in 
close analogy to Eq. (pSj) we then find that the modified kernel Kk{x;T]) can be 
rewritten as 

Kk{x;r]) = Kkix;r])-Fkix;7]), (44) 

where for even and odd k, respectively, r]''^'^Kk{x;r]) is given by Eqs. (IATBD and 
( [A 15| ), if in the integrals Zk is replaced with Zk{ri) =Zk/rj. Finally, Eqs. ( pQ]) and 
(|30D , respectively, are replaced with 

2a;(m-|-l)! 



2m+l 



'x;ri) 



(27r/?7)™+3/2 7o \ ^ ' 

r2™$[m+2, 3/2, -x'^\{r'^; r]) tanh(rV2)] 
e-(rn+i)r2 sinh'"(rV2)[A(r2; r]) cosh(rV2)] 



m+2 



(45) 



and 



Kim{x\r]) 



X 



' ' (27r/r/)"^+i 
„2m-igmrV2 1/2, -x^Xir"^] T]) tanh(rV2)] 



sinh™(r2/2) L e-('"+i)^'-V2[A(r2;?7)cosh(r2/2)]"+i 



where 



X{r^;ri) = 1 + (r/ - 1)(1 + e"'-')" 
and the polynomial Fk{x; rf) reads as 

~'-2r]Y{k - n) 



(46) 



(47) 



Ffc(x;r/) 



E 



oo 
1=0 



Hfc-2n(a;) 
1 



(Z + 1) ...(/ + A;) 



(48) 
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Needless to say that the polynomial can again be omitted since both Kk{x;ri) 
and Kk{x\ rj) are solutions of the problem and in practical measurements Kk{x; 77) 
is more suited for error reduction than Kk{x;ri). The numerical calculation of 
Kk{x; rj) can be performed in a way as outlined in Sec. ^ for Kk{x). Examples of 
Kk{x] rj) for various values of k and r] are shown in Fig. 0. 

It should be pointed out that the sum rules ([A 13| ) and ( |A 14|) used in the 



derivation only apply when |2fc(^)| < 1- Hence we observe that the condition 
r) > 1/2 must be fulfilled in order to compensate for Gaussian data smearing, 
which is analogous to the density matrix reconstruction in the Fock basis It 
is worth noting that the condition > 1/2 corresponds to the requirement that 
the width of the Gaussian (^) is smaller than the vacuum noise. From Fig. ^ 
we see that with increasing |x| the kernel Kk{x\rj) for odd k rapidly approaches 
the classical limit ( p!3|) for perfect detection, whereas for even k it approaches the 
classical limit (0) up to an irrelevant ?7-dependent constant. The results reveal 
that in classical optics it is impossible to compensate for the losses in nonperfect 
detection, because of the vanishing vacuum noise of a classical oscillator. 

As expected, substantial differences between Kk{x; rj) and Kk{x) are observed 
in the region around a; = 0, and they increase with decreasing r] (Fig. |^). The [com- 
pared with Kk{x)] stronger variation of Kk{x]rj) implies that the use of Kk{x;rj) 
for sampling of the exponential phase moments from the smeared quadrature- 
component distribution gives rise to a larger statistical error than the use of 
Kk{x) (for the statistical error, see Sec. ^I3|) . This is obviously the price paid 
for suppression of the systematic error. Based on the precision of the data avail- 
able, the experimenter should therefore decide whether to use Kk{x;r]) (which 
increases the statistical error) or K^^x) (which decreases the statistical error but 
introduces a bias). 

4.2 Phase discretization 

In practice, p{x,{}) can only be measured at A^ discrete phases "i?;. When they 
are equidistantly distributed over a 27i interval, i.e., {}i = {2n/N)l, where / = 
0, 1, . . . A^ — 1, then application of Eq. (^ yields the experimentally determined 
exponential phase moments 

^fc(^) = TfT. e'""' / dxKk{x)p{x,^i), (49) 



which can be rewritten as, on using Eq. (|23|), 

N-l 00 



^fc(^) = E E e'-^'-^^^'g^,n / dxKk{x)^Ux)Mx)- (50) 



1=0 m,n=0 
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Taking into account that A^~^ E/Io^ e*2^(*^+"-'")'/^ = 4+n-mmodiv and recalling 
Eqs. (|0D and (|2|), we derive 



^k{N) = ^k + A^'^^\ (51) 



where 



oo oo 



~ 51 51 {(^ri+k+sN,nQnlk+sN,n + Qn,n+sN-kQn}i+sN-k) (^2) 



=ln=0 



represents the systematic error owing to phase discretization. In Eq. ( ]5^ the 
abbreviation 

/■oo 

Q)^1 = 2tt dxKk{x)i)m{x)i)n{x) (53) 

is used and it is assumed, for notational convenience, that N >k. Note that from 
physical arguments it is also reasonable to assume that the number of phases is 
larger than the index of the measured moment (otherwise the systematic error 
could dominate the result). 

From Eq. (^2]) we see that the error is influenced by all off-diagonal density- 
matrix elements of the type of Qn+k±sN,n- The effect, which is also called "alias- 
ing" , has also been found in the reconstruction of the density matrix in the Fock 



basis from the data measured in balanced [|19[ and unbalanced homodyning. 
For highly excited states (i.e., ^n,m~0 if n,m<no, with rio ^ 1) the relevant Qmn 
can be approximately calculated, using in Eq. (^) the classical kernel given in 
Eqs. (|TD and (0): 



n(^) ^( ^ o^^^ ^ -O^'') r54l 

^n+k+sN,n \ ^) gJ\f_^J^ ' ^n,n+sN-k ^n-k+sN,ni \^^) 

Note that sN is even, since from the symmetry properties of Q^^nn follows that 
Qnin = Oifm + r;, + /cis odd. Combining Eqs. (|5^) and (^4]) yields 



for odd and 

A*f.|,-l)-/^(-J^*..„. + -JL_,_) (56) 

for even N . We see that the error of the measured exponential phase moment is 
expressed in terms of higher-order moments, and it decreases with increasing N . 
Note that the difference between the errors in Eqs. ( ^5|) and (^) reflects the fact 
that with regard to a vr interval, the number of different phases is twice as large 
for odd as for even N , which of course substantially reduces the systematic 
error in the first case. 
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4.3 Statistical error 



When in an experiment n{i}i) measurements are performed for each phase -t^i, the 
exponential phase moments can be estimated as, on applying Eq. (||), 



2tt 



N-l 



N 



1=0 



J2 Kk[xMl 



I) r=l 



(57) 



where Xr{'&i) is the result of the rth individual measurement at the phase -^i. 
Taking the average of all estimates ^l''*^(iV), 



(=0 



271 



N-l 



JM, 



yields, as expected, "^k^N) from Eq. (|i9|). 



(5^ 



k ! 



because of 



KkiXrm) 



dxKk{x)p{x,di). 



(59) 



(60) 



The variances of the real and imaginary parts of '^^I^^^\n) can be obtained in a 
similar way. Taking into account that the individual measurements are indepen- 
dent of each other, we derive 



and 



Var{Re[^f*^(Ar)]} 



Var{lm[^f*^(Ar)]} 



47r^^^cos^ 



'-^ {/\Kk[xr{di)]y 



(61) 



where 



{^Kk[xrmy 



dx KI{x) p{x,'&i] 



dx Kk{x) p{x,{}i 



(62) 



(63) 



Equations (pT| ) - ( p3| ) enable us to estimate the statistical error of the mea- 
sured moments, substituting in Eqs. ( [6T| ) and ( p^ ) for {AKk[Xr{'&i)]}'^ the cor- 



responding estimates. From Eqs. ( pTD - (63) it is seen that the statistical error 
depends on the shape of the function Kk{x). In order to reduce the statistical 
error, the ambiguity in the determination of Kk{x) can be advantageously used 
to choose it such that it varies as slowly as possible. This is one of the reasons 
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for omitting the polynomial in Eq. (pHI). Moreover, the statistical error can be 
reduced when the number of events, n{{}i), is appropriately varied with the phase 
From Eqs. (BTI) and (ET) it is suggested to increase n{^i) for such phases for 



which {AKk[xr{i^i)]}'^ , Eq. (|63|), becomes relatively large. This is typically the 
case when p{x, i)) is essentially nonzero in an x interval around a; = 0, in which 
Kk{x) strongly varies with x. Note that this result is in qualitative agreement 



with that of the maximum-likelihood method for estimating phase shifts [21 



5 Computer simulations of measurements 
5.1 Exponential phase moments 

To illustrate the method, we have performed computer simulations of measure- 
ments of the quadrature component distribution p{x,{}), assuming the signal 
field to be prepared in various states, such as a squeezed vacuum |0)s = 5'(^)|0) = 
exp{ — ^[^{d^)'^ — ^*a^]}\0) and a displaced Fock state \a, n) =D{a)\n) =exp{ad^ — 
a*d)\n). We have restricted attention to perfect detection and assumed that the 
measurements are performed at A^ = 120 (equidistant) phases -t^i within a 27r inter- 
val and n{'di) = 10'^ events are recorded at each phase. Examples of the sampled 
exponential phase moments are shown in Figs. ^ and ^ for a squeezed vacuum 
and a displaced Fock state, respectively. The error bars indicate the standard 
deviations obtained according to Eqs. ( |6TD and (^). Compared to the statistical 
error, the systematic error due to phase discretization is (for the chosen number 
of phases ^i) negligible small. From Figs. ^ and ^ we see that [for the chosen 
numbers n{i)i) of events] the exponential phase moments are obtained with suffi- 
ciently good accuracy, provided that k is small enough. We further see that the 
accuracy decreases with increasing k. (Note that for the chosen state parameters 
the imaginary parts must vanish for all k.) Clearly, the accuracy can be improved 
if the number of measurements is increased. 



5.2 Phase distribution 

The possibility of direct sampling of exponential phase moments offers novel 
possibilities of experimental verification of fundamental number-phase uncer- 
tainty relations, as has been shown recently [^. It is worth noting that the 



measurements can be performed with high precision since only low-order mo- 
ments play a role. Here we address the problem of the determination of the 
whole phase distribution P{(p). 

Since the exponential phase moments are nothing but the Fourier compo- 
nents of the phase distribution [Eq. (^], the sampled moments can be used to 
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reconstruct the phase distribution according to 



1 

2^ 



(64) 



Moreover, since any physical quantum state can be approximated to any desired 
degree of accuracy by truncating it at some photon number n^ax if ^max is suitably 
large, from Eq. (^) it follows that (for chosen accuracy) the number of moments 
"^k that effectively contribute to P{(p) in Eq. (|6^ ) is finite, i.e., \k\ = l,2, . . . , K, 
with K = n^s,x- Hence, P{(p) can be obtained truncating the sum in Eq.(^) at 
\k\ = K and substituting for the the measured moments ^'if''*''. The phase 
distributions that are reconstructed from the measured moments given in Figs. ^ 
and I for a squeezed vacuum and a displaced Fock state, respectively, are plotted 
in Fig. 1^, on assuming that nmax = 20. The statistical error of ^1/^^^*^ gives of course 
rise to an error of P{^p)- Since the error in P{ip) can be obtained easily from the 
law of error propagation in a standard way, we renounce the calculation here. 

Finally, it should be pointed out that there are other methods, such as least- 
squares inversion |2^, and maximum-entropy inversion |]2^ , which can be used 
for reconstructing the phase distribution from the measured (i.e., inaccurate) 
exponential phase moments - methods that have been successfully applied in 
various fields of physics. Let us briefly comment on the application of the method 
of least-squares inversion. For this purpose we return to Eq. (^ and ask for P{if) 
that best fits the experimental data at M chosen phases (m = 0, 1, . . . , M— 1, 
with M ^ 2K). An answer can be given applying the method of least-squares 
inversion |2^, ^ to the set of 2K linear equations for M unknown P{ipm) {<^m = 
27rm/M), 



2-K ^"^-^ 

m=0 



27r 



A/-1 



(A; = 1, 2, ... , K), i.e., minimizing the functional 



(65) 
(66) 



X 



k=l 



(Re) 
k 



2n 



M-l 



+ 



a 



(Im) 



27T *^ 



M-l 



Im - V7 II P{^m) sm{kipr. 

-'^^ m=0 



(67) 



where a'^^'' and a 



(Im) 



respectively, represent the errors involved in the determi- 



nation of Re "^^^^^^ and Im \E'^'^^*'' . In particular, when cx^^^^ ~ o"!^™"* = then (for 
M^2K) the resulting P{ifm) is in agreement with that obtained from Eq. 
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with I A; I <K. From comparison with the exact phase distribution it can be seen 
(Fig. H) that outside the regions in which the phase distribution is essentially 
nonzero artificial oscillations and even negative values of the reconstructed dis- 
tribution may be found, because of the statistical error of the measured moments 
(cf. Figs. ^ and The problem can be partially overcome introducing regular- 
ization techniques in the solution of Eqs. ( |65D and (|66D (for details, see p3|), as it 
can be seen from Fig. The figure also reveals that the artifacts are suppressed 
at the expense of a smearing of the overall distribution. 



6 Summary and conclusions 

We have presented a method for direct sampling of the exponential moments 
of the canonical phase of a single-mode radiation field from the data recorded 
in balanced homodyning. The sampling method enables us to determine the 
moments in real time, together with the statistical error. It is worth noting 
that the method renders it possible to closely relate the basic-theoretical con- 
cept of canonical phase to the experiment. The sampling functions relating the 
quadrature-component distribution to the exponential phase moments are well 
behaved. With increasing quadrature-component they rapidly approach their 
classical counterparts given either by step functions (for odd moments) or loga- 
rithmic functions (for even moments). In this way the concept provides us with 
a unified approach to the experimental determination of the canonical phase in 
both quantum and classical optics. 

In our approach to the construction of the kernel functions needed for direct 
sampling the exponential phase moments we have extended the proposal made 
in [jl8l . Hence the here used subtraction of the polynomial arbitrariness from the 
sampling functions can be considered as a significant improvement. In fact the 
omittion of the ambiguity from the kernels provided us with functions that can 
be universally used for any physical state and that are much more insensitive to 
errors. This point should be also considered if a similar approach to the derivation 
of the kernel corresponding to any other generic field operator is adopted. 

In order to study the accuracy of the method, we have discussed the infiu- 
ence of various experimental inaccuracies on the measured exponential phase 
moments. In particular, the finite number of local-oscillator phases results in 
an aliasing effect, whereas smearing of the quadrature-component causes a bias 
toward smaller absolute values of the moments. When the data smearing results 
from imperfect detection with efficiency larger 50% 77, then modified sampling 
functions can be introduced for compensating the losses. 

To illustrate the applicability of the method, we have performed computer 
simulations of measurements for two quantum states and and presented the mea- 
sured exponential phase moments including the statistical error. Finally, we have 
used the moments for a reconstruction of the whole phase distribution. The re- 
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suits obtained are in good agreement with the theoretical predictions. 
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Appendix A Derivation of Eqs. - (CT) 

Combining Eqs. (]T9|), (|25|), and ( p6D yields, provided that the moments and cor- 
relations (a^"a™) exist for all n and m, 

^k= dxKk{x)e''"^p{x,^), (Al) 

J2-K J —oo 



where 



oo 



K,{x) = {27T)-'J2cr^2i+k{x), (A 2) 

1=0 

with 

Using the relation 

1 1 f + CO , f + (X> , „ „ 

dtie-*? ■■■ / rftfcc-'^*^ e-"'-s (A 4) 



'{n + l)...{n + k) vr'^/^J-oo 

where 



k 

2 \ " ,2 



E^?' (A 5) 



we may rewrite Eq. ( |A 3|) as 



(27r)'"+5(2m + 1 + 2/)! ^-oo 

c^Wie-(2'"+^)*-+^zL+i (A 6) 



and 



(2^) ^ (2m + 0! /•+ 
^ ~ (27r)'"(2m + 2/)! J- 



rftie-*? ■■■ / rfWe-2"^*^'"4^ (A 7) 

-oo J —oo 
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for k = 2m+l and k = 2m, respectively, where 

z, = \{e-< - 1). 



(A 



We now substitute in Eq. (Q) for the expressions (|A 6|) and ( [A 7] ) and 
change the summation index / as m + / = j. In order to separate from Kk{x) an 
irrelevant polynomial Fk{x), we decompose Kk{x) into two parts, 

Kk{x) = Kk{x)-Fk{x), (A 9) 

on rewriting the j sum such that in Fk{x) it runs from j = and j = l, respectively, 
to j = m — 1 for odd and even k. In this way we find that 

(m + 1)! 



2m+l 



X] 



(27r)"^+3/2 



X 



-j-oo g 
dt; 

oo 



dti e 



-(2m+l)t| 



and 



(2vr; 



m+l 



dti e~ 



X 



+00 g~2mi^^ 
dt2m 

-oo 



^2m 



E 



r(m + l+j>^„ 



(2j)!r(m+l) ^ 

The polynomial Fk{x) can be written as 
1 



(All) 



(27r)i+^'/2 



+00 g-*^*! 
-oo Zj, 



(A 12) 



where the summation index j has been changed as m — j = n. Next we apply the 
sum rules ||2^ 

1 



oo 



°° Tin -I- iW-? 



2x 



-<l>[a,3/2,xV(l + ^)] 



(A 14) 



^t^o(2j + l)!r(a) ^^^^ ^ (1 + z)^ 

[\z\ < 1; $(a, 6, y), confluent hypergeometric function] to the j sums in Eqs. ( [A IQ ) 
and ( [A 11|) and obtain 

2x(m + 1)! /■+°° 



-f^2m-f 1 (a;) 



(27r)™+3/2 



(iti e 



X 



+ 00 



dt 



2m+l e 



'(2™+l)tL+i 



^[m+2, 3/2, 2:2>n+i(l + ^2m+i ) ^a;^^ 



KA 15) 
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hoc 



J^2m{x) = -, r -r / rfti 6 ' 

^ (27r)'"+l i-oo 



X . . . r e-L { '^^'']!^r^^:^^^ - 4r} (A 16) 

To write the multidimensional integrals in Eqs. ( |A 15|) and ([A 16| ) in a more 
compact form, we change the variables as, on using generalized spherical coordi- 
nates 

i-\ k-l 

ti = r cos (pi sin ipj if i < k, and tfc = r sin (pj, (A 17) 

3=1 3=1 

and r = rk, with <r < oo, 0<ipj <tt if j < A;— 1, and < ipk-i < 27r. Introducing 
the function 

Q{k) (^2) _ r g-t? r d^. e-^*? sin'^-^-i 

Jo Jo 

dipk-ie~^^-^^''-^e-'''l (A 18) 



X 



2n 







and recalling Eq. ( |A 8| ) [together with Eq. ( |A 5| )], we find that Eqs. ( A 15| ) 



and ( |A 16|) can be rewritten as given in Eqs. (129|) and (0) in Sec. |^. Finally, we 
expand in Eq. ( |A 12| ) 



and perform the Gaussian integrals to obtain Ffc(x) in the form of Eq. ( |3TD in 
Sec. S. 



Appendix B Derivation of Eqs. and ( [S3l ) 

In order to write the function il^'^^z) in the form of Eq. (|32D together with 
Eq. (|33D, we first rewrite Eq. ( |A 18| ) as 

Q{k) (^^2^ ^ r r sin''-'-^ ipi--- r e-^" , (B 1) 

Jo Jo Jo 

where 

Xk=tl + 2tl + --- + ktl, (B 2) 



tk being given in Eq. ([A 17| ). We then expand the exponential e "^'^ in a power 



series, which implies the power-series expansion of 

oo 



m=0 
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Here, the expansion coefficients are given by 



d(pi sin'^ ^ (pi 



d^, sin^-^^Vi-- - / d^k-iY,"', (B4) 
Jo 



2n 



m\ Jo 

where Yk = Xk/r"^ is an r-independent angular function, Yk= Yk{(pj). Particular 
integrals of the type given in Eq. (P 4|) are calculated in They can be used 
to prove, by induction, that for arbitrary k and m the coefficients Ajf^^ can be 



written in the form given in Eq. (33). 



Appendix C Asymptotics of Kk{x) 

In order to find the asymptotic behaviour of Kk{x) for large x, we change the 
variables in the integrals in Eqs. (^) and ( |30D according to tanh(r^/2) =7 and 
obtain 

K2m+i{x)=x ['d^D(^"^'\^)^m + 2,l,-jx^) (C 1) 

Jo 



and 



where 



K2mix) = d^D^^^\^) [$(m + l,i,-7a;2)-(l+7)- 



(m+l) 



2(-l)"^(m + 1)! 



(27r)'»+3/2 



In 



In 



1+7 
1-7. 



1+7 
(I + 7) 



-1 m-1/2 



2m+l 



7^(1 - 7) 



(C 2) 



(C 3) 



D(-)(7) = 

^ ^' (27r)™+i 



In 



m— 1 



X 



In 



1 + 7' 
1-7. 



1+7 

,1-7. 

(I+7) 
7™(1 -7) 



2m 



(C4) 



Let us ffist draw attention to K2m+i{,x) in Eq. ( |C 1|) . For 7 > the function 
^{2m+i)j^^^ is finite, and from the asymptotic behaviour of Q^-^^r"^) for large r 
it follows that D^'^'^^\l) = 0. From the asymptotic behaviour of the confluent 
hypergeometric function. 



$(a, c, -t) 



r(c) 1 

r(c-a) t"^ 



t — »• 00, 



(C 5) 
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it follows that $(m + 2, |, — 7x^) — if |x| — oo except for 7 = 0. Hence, for 



^ 1 we can approximate the integral in Eq. (|C 1|) , replacing D^'^"^^\'y) with 



its expansion for 7<^ 1. Taking into account that 



In I 



1 + 7 
1-7. 



27 1 + 0(7' 



(C 6) 



and recalling the expansion of n^''\z), Eqs. (§3) and (|B|), we find that 



In 



1 + 7 
1-7. 



T{k/2) 



l-(A; + l)7 + 0(72) 



Thus, for 7<^ 1 the function D^'^"''^^\'j) , Eq. ( |(J 3| ), can be given by 



(-l)'"(m + l)! 
7rr(m+ 1/2) 



7 



^1/2 



1 + 0(7') , 7«1, 



(C 7) 



(C 



so that Eq. ( |U l\ j for x ^ 1 can be written as, on changing the variables as j^^'^x 



2(-l)™(m + l) 



dt 



l + 0{tyx^)] <l>(m + 2,f,-t' 



7rr(m + l/2) Jo 



7rr(m+ 1/2) 
where the abbreviating notations 

/(i)(x) = / rft<l>(m + 2,f,-t2) 

J X 

and 

I^^\x) = / rftt^$(m+2,|,-t2) 



(C 10) 



(Cll) 



have been introduced. The integral I^^\0) can be calculated to be ||26 

7r(2m + 1 



/«(0) 



rft<l>(m + 2, |,-t2) 



fm + 1)! 



r[(2m+l)/2] (C 12) 



and it can be shown that /(2)(0) =0. Further, from Eq. ( FSj) it follows that I^^^x) 
= 0(x-''"+^) and /('^(x) = 0(x-'™+^). For x < (|x| > 1) the calculations are 
quite analogous, so that 



K2m+i{x) = i(-l)™(2m + 1) sign (x) [1 + 0(x- 



2m-3\ 



, |x|>l. (C 13) 



To find the asymptotic behaviour of K2m{x), we subdivide the interval of 
integration in Eq. ( |C 2|) as 



K2mix) = J(0, M/x') + J(M/x', 70) + /(70, 1), 



(C 14) 
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where 



I{a,b)= f d-fD^^"'\-f) [<l>(m + l,i,-7x2)-(l+7)— 1 



(C15) 

and 7o <^ 1 and M ^ 1 such that M/x^ < 70. This reflects the quahtatively 
different behaviour of the integrand in these intervals. In the integrals /(O, M/x^) 
and liMjx^ ^ 70) the variable 7 is small, so that L>(2m)(^^ 

can be expanded as, on 

using Eqs. (JCJ) and (|CTD , 



D(2-)(^) = (27r)-i(-l)"^m7-i [l + 0(7')] , 7 < 1. 



(C 16) 



We also expand (I + 7) ^ and can rewrite J(0,M/x^) as (after changing the 
variables) 



/(O, Mjx^) = Ci{M) + 0(m/x^) , 



where 



Ci(M) = (2 



:7r 



M dt r 

T 



^(m+l, |,-t) -1 



(C17) 



(C 18) 



J(M/x^, 7o) and /(70, 1) can be calculated integrating the two terms in Eq. (|C| 
14| ) separately. Recalling Eq. ( P 5| ) and changing the variables, we find that 
/(M/x^, 70) = + J(^), with 



= (27r)-i(-l)'"m / [1 + 0(tVa;')] $(m + l, i, -t) 



= o 



(70a; 



'A/ t 
+ C2(M) 



(C 19) 



(27r)-i(-l)'^m 



1 + 0(f )] [1 + 0(7)] 



where C2(M) = C(M"'"-i), and [after expanding (1+7)-"""^; 

no (i7 

/Af/x2 7 

= 7r-i(-l)'"+^m In |a;| + C^iM) + ^4(70) + 0(m/x' 
where 

C3(M) = 7r-^(-l)'"+imlnM, 

^4(70) = (27r)-i(-l)'"+im [In 70 + O(7o)] • 

Finally, /(70, 1) can be written as, on using Eq. (|C 5|) , 

J(7o,l) = /'rf7^^''"n7){(^[(7^')-"-']-(l + 7)"™"'} 
= C,{^o) + O{%"^x-'"^-'' 



(C 20) 



(C21) 
(C 22) 



(C 23) 
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where 



C5(7o 



d^Z}(2™)(^)(l + ^)- 



m—1 



(C 24) 



70 



Substituting in Eq. ( |C 14| ) for I{0,M/x^), I{M/x'^,'jo) and /(70, 1) the expres- 
sions derived above, it can be shown that the logarithmic divergences of 6*4(70) 
and 6*5(70) cancel for 70 —>0 and those of C3(M) and Ci(M) cancel for M-^ 00, 
and hence 



K2m{x) = 7r-i(-l)™+im In |a;| + 0[m/x' 



+0 (7ox' 



-m— 1 



C(7o,M), |x|>l. 



(C 25) 



Here, the constant C(7o, M) is given by the sum of the constants Ci - C5. Note 
that for 7o ^ and M — > 00 the constant C(7o, M) becomes independent on 70 
and M [for a determination of this irrelevant constant the value of the integral 
(|C 24|) must be known]. 
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Figure 1: The x-dependent part Kk{x) of the samphng function Kk{x^d) = 
Q^k^Xk{x) for the determination of the exponential phase moments from the 
quadrature-component distribution i?) is shown for various odd (a) and even 
(b) k. 
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Figure 2: The function gk{x;ri) for the determination of the systematic error 
A*^'^)\l/fc, Eq. (^OD, which is associated with Gaussian data smearing, is shown for 
k = l,2 and rj = 0.6. 
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Figure 3: The x-dependent part Kk{x;r]) of the samphng function Kk{x,'d]r]) 
— e^'''^Kk{x;r)) for the determination of the exponential phase moments from 
the smeared quadrature-component distribution p{x, i); rj) is shown for various k 
and rj [r] — l (full lines), 7^ = 0.8 (dashed lines), 77 = 0.6 (dashed-dotted lines)]. 
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Figure 4: Examples of measured exponential phase moments ^fc'^^*'' are shown 
for a mode prepared in a squeezed vacuum |0)s = S'(^)|0), with ^ = —1.31, i.e., 
(n) = 3 [(a) real part of (b) imaginary part of ^fe^***'']. The error bars 

indicate the estimated statistical error. In the computer simulation it is assumed 
that the quadrature component distribution p{x, -§) is detected at 120 phases ■»? 
equidistantly distributed in a 27r interval and that at each phase 10^ events are 
recorded. 
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Figure 5: Examples of measured exponential phase moments "^^^^^^ are shown 
for a mode prepared in a displaced Fock state \a,n) = D{a)\n), with a = —1.5 
and n = 2, i.e., (n) =4.25 [(a) real part of ^[f'**'*; (b) imaginary part of ^'^'^***'']. The 
error bars indicate the estimated statistical error. In the computer simulation it 
is assumed that the quadrature component distribution p{x, -§) is detected at 120 
phases ■»? equidistantly distributed in a 27r interval and that at each phase 10^ 
events are recorded. 
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Figure 6: The canonical phase distribution P{<f) reconstructed from K = 20 
measured exponential phase moments "^^^^^ given in Figs. ^ and |^ is shown for 
the squeezed vacuum (a) and the displaced Fock state (b) therein. The results of 
direct application of Eq. (p^ ) (dashed-dotted lines) and application of regularized 
least-squares inversion (dashed lines) are compared with the exact distributions 
(solid lines). 
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